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Abstract 

A model of quantum noisy channel with input encoding by a classical 
random vector is described. An equation of optimality is derived to de- 
termine a complete set of wave functions describing quantum decodings 
based on quasi-measurements maximizing the classical amount of trans- 
mitted information. A solution of this equation is found for the Gaussian 
multimode case with input Gaussian distribution. It is described by the 
overcomplete family of coherent vectors describing an optimal quasimea- 
surement of the canonical annihilation amplitudes in the output Hilbert 
space. It is found that the optimal decoding in this case realizes the maxi- 
mum amount I = Sp In [l + Sy(iV-|-l)] of the classical information as trans- 
mitted via the classical Gaussian channel with the effective noise covari- 
ance matrix N + I. A physical realization of optimal quasi-measurement 
based on an indirect (heterodyne) observation of the canonical operators 
is suggested. 



Introduction 



The use of homodyne or phase-sensitive receivers for the reception of a coher- 
ently modulated quantum signal enables one to extract in the quaiclassical limit 
only a half of the encoded information. The coherent reception of quantum sig- 
nals, based on simultaneous direct measurement of amplitude and phase (or 
phase coordinates of each mode of the signal), is forbidden due to the uncer- 
tainty relations of quantum theory. Therefore, it is of interest to investigate the 
potential possibilities of decoding in quantum communication channels based 
on the indirect, heterodyne methods of observation. The theory of indirect 

'Originally published in: Radio Eng. Electron. Phys., 19 (9), p. 1349, 1973. [trans, from 
Radiotekhnika i Electronika, 1973, 19, 9, 1839. 844] 
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measurements of incompatible observables {bj} has been under the develop- 
ment starting from the pioneering studies 0J [7] to meet the demands of the 
quantum communication theory. Using a quantum reformulation of the sta- 
tistical inference theory for a Bayesian risk criteria it was recently shown how 
to formulate optimal quantum information processings based on the indirect 
measurements in the classical-quantum noisy channels. 

The use of the concept of indirect measurement has made it possible to 
indicate for a Gaussian classical-quantum communication channel the amount 
of information 




transmitted by each mode of coherently modulated signal in using a quantum 
receiver with a large linear amplifier* (s, n are the energies of the corresponding 
modes of the classical input signal and of the quantum additive noise expressed 
in units of his). The same amount of information was deduced in [S], where a 
method of its decoding was indicated, leading to indirect measurement. The 
rate of information transmission by the coherent signal was also determined in 
an earlier work [B], where vectors of the {'fig} coherent states were used as the 
operators characterizing the coherent reception. The problem of the realization 
of the coherent measurement and its optimality was not considered. However 
the natural question - whether the decoding based on the quasi-measurement 
described by the coherent vectors (|1.4(l is optimal, and if so, in what conditions 
it is unique, has not been resolved. 

In order to answer this question we derive the equation determining the op- 
timal decoding vectors {'fig} according to the criterion of the maximum amount 
of information and specify the optimal indirect measurement realizing the cor- 
responding quasi-measurement. 



1 Decodings based on quasi- measurements 

Quasi-measurement of noncommuting operators {bj} with values in a given 
(measurable) space B 9 (5 are described by operator-valued (positive, cx-additive) 
measures n(d/3) normalized to the identity operator 1 in the Hilbert space H of 
quantum states as II (B) — 1 such that 

i = J n(d/3), J p 3 u(dp) = b r (i.i) 

If the resolution of identity in (|1 . 1|) is orthogonal such that II(d/3)II(d/3 ) = 
II(d/3 n d/3 ), then it describes the usual (direct) measurement of compatible 
observables in the spectral decomposition bj = f f3jH(d/3). In the general case, 
however, the operators {bj} do not have common orthogonal spectral decompo- 
sition, and it can be said that the nonorthogonal expansion (|1.1|) describes their 

"This result was presented by R.L. Stratonovich at the symposium on information theory 
(Dubna, June 15-25, 1969). 
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simultaneous quasi-measurement. The probabilistic operator measure H(d/3) is 
an analog of classical randomized decision rules P(df3\b). 

The classical probability distribution P(d/3) of the results of quasi-measurements 
in the state specified by the statistical density operator p is determined by the 
Hilbert space trace Tr of its product with the quantum probability measure 
II(d/3) in the same way as in the case of direct measurements: 

P(df3) = Trn(d/3)p. (1.2) 

Below we shall specify the operator II(d/3) in the form 

n(d/3) = <p p <p}n(d0), (1.3) 

explicitly taking account of its positive defmiteness and cr-additiveness. Here 
{(fp,/3 £ B} is a family of non-Hcrmitian operators acting in general from a 
Hilbert space F into the original space H and, together with a positive mea- 
sure /x(d/3) as the reference measure on the set B, satisfying the condition of 
normalization (| 1 - If) as the completeness condition for {</?g}. In particular, if 
{ifp} are operators acting from one-dimensional complex space F = C into H, 
i.e., are a set of wave functions (fa £ H (in general nonorthogonal and non- 
normalized: tpgippi 7^ <^8/3')j t nen the quasi-measurement specified by them is 
called elementary. We note that any quasi-measurement can be realized by 
direct measurement of a certain family of compatible observables ^Pjj m an 
expanded quantum system consisting of the original system and an independent 
auxiliary system. This follows from the well-known fact that any operator mea- 
sure n(d/3) can be represented as partial averaging H(d/3) = Trn 6(d/3)/9 of an 
orthoprojector measure 0(d/3) given by common eigenvectors of the commuting 
operators $j = J (3jQ(df3) acting in the expanded space H £g) Hd, along with 
the state p of the auxiliary system described by a Hilbert space Ho. Such a 
realization of quasi-measurement is called indirect measurement. 

As an example we consider an elementary quasi-measurement specified by 
the coherent vectors 

<p a =exp | -^a+a + a+aj |0) = |a). (1.4) 

It can be realized by direct simultaneous measurement of the complex ampli- 
tudes 

& = a + a£ = (&j), j = l,...,r, (1.5) 

where a — (otj) and a = (aj) are columns of the amplitudes aj G C and the 
annihilation operators aj operating into the original Hilbert space H, ot* and a^ 
are the adjoint raws, while a$ denotes the column (a^) of creation operators 
a^- of an auxiliary (heterodyning) system in the vacuum state 1 0) o £ Ho. For 
this purpose it is sufficient to note that the partial averaging over the vacuum 
state |0)o amounts to the determination of the projections <p a of the eigenvectors 
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tp a € H ® Ho of commuting operators a = a + Oq in the original space H. One 
can easily verify that the generalized vectors 

?A Q =exp{-( a -a)ta+}|0)o«)|a} (1.6) 

satisfy the equations aip a = aip a , a*ip a = °*^a an d f° rm a complete orthonor- 
mal set with the measure 

r 

d/z(a) = Y\ ^dReajdlmaj. 

3=1 

Multiplying these vectors from the left by o(0| we find that the projections 
<p a G H are actually the coherent vectors \a) and, together with the measure 
/i(da) = d/i(a), determine ideal measurement corresponding to the indirect 
measurement of noncommuting aj. 

The question - whether the decoding based on the quasi-measurement de- 
scribed by the coherent vectors (|1.4|) is optimum and if so in what conditions 
- is of interest. In order to answer this question we are going to derive the 
equation determining the operators tpg, which are optimal according to the cri- 
terion of the maximum amount of mutual information and specify the optimum 
quasi-measurement. 

We shall assume that the family of statistical operators |p(??)} defining the 
state of the quantum communication channel is given as a function of the random 
information parameters $ = having the distribution P(dz9). The amount of 
Shannon information of the parameters d and the results of quasi-measurements 
(3 is determined by the usual formula 

!/*,* = jj pW)ln^^P(dtfMd/3), 

where according to Q1.2|) . the density p(f}\&) has the form 

= Tr v^ppifi) = Tr ^p(tf)^, (1.7) 
p{0) = / p P(di9), satisfying the normalization condition 

j p(p\0)fi(d0) = l, J pfflridp) = l. 

Using the method of Lagrangian multipliers we set up a function 
J W)Tr^p(tf)^P(dtf)-Tr^A^J ^(d/3), 

where i((3,d) — In [p(/3|i?) / / f>(/3|$)P(d#)] , and A is an operator determined 
from the condition of completeness J ipp(p%/j,(d/3) = 1. Varying it over tp% we 
obtain the equation determining optimum tpp. 

(j(/3)-A)^=0, (1.8) 
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where I (f3) = J 0)p(0)P(d<?) and A = J I ((3) <p <p* M (d/3) . 

In the derivation of this equation we have made use of the fact that 

J P (/3\tiMMM#)Pm = J p(W)6v* (jrjPL _ -P-\ ^P(dtf) = 0, 

where p — J* /?($)P(d$). Multiplying equation l|1.8(l from the right by <p% and 
integrating with the measure ^(d/3) we determine the operator A: 




P 0?)i(M^/x(d/?)p(dtf), 



the trace of which Tr A gives the maximum amount of information. And, finally, 
eliminating the operator from Equation. (|1.8|) and multiplying it from the left by 
A, we rewrite this equation in the following equivalent (because of the condition 
of completeness J <p a tp* a p.(da) = 1) form: 

j <p{$) U m *)- J <Pp<p%.Vfi\(P , P{M) = 0. (1.9) 

The equation thus obtained is a complex nonlinear equation in (fg and it is not 
possible to solve it explicitly in the general case. However, as will be shown in 
the next section, in the Gaussian case the coherent vectors ()1.4J) are optimum 
operators <pg. 



2 Optimal decoding in quantum Gaussian chan- 
nel 

Let the received signal b be a superposition b = d + a of a complex vector 
= = 1, . . • , r) and an independent Gaussian Boson noise a — (aj) (aj 

are annihilation operators of the investigated modes: [aj, — Sjk, having zero 
mathematical expectation (aj) — and specified average number of quanta 
(a'j'aj) — fij). In Glauber's representation [2] the statistical operator of this 
noise has the form 

p = J \a){a\\N\- 1 exp{-a' i N- 1 a}dp(a), (2.1) 

where \a) is the coherent vectors (II .4|) : N is a matrix whose eigenvalues are rij, 
and |JV| = det N = n\ ■ ■ ■ n r . In formula (|2.1(l the matrix form of writing scalar 
products is used, a' r N~ 1 a = J2j k a *j(N~~ 1 )jkUk , similar as a^a = J2j a t a i 
in 1)1.4(1 . The signal b after passing through the indicated linear channel is 
described by a family of operators {/?($)} of the form 

p(0)= / \a)(a\p(a\d)dp(a), (2.2) 
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where p(ai|i9) = | TV^ | 1 exp { — (a — , &)^N~ 1 (a — $)}. The probability distribution 
P(di?) of the transmitted signal chosen from the condition of maximum en- 
tropy and finiteness of the energy S~ 1, d) < (S is the given correlation matrix 
(flj'dk)), also has a Gaussian form: 

r 

P(dtf) = |S'|- 1 exp{-i9 t S'- 1 ^}d M (i?), d/z(i?) = JJidRe^dlm^. (2.3) 

In the case under investigation it is natural to expect that the results of 
optimum quasi-measurement of the Gaussian observable b also have a Gaussian 
distribution; among the Gaussian operators (p^, determining this distribution, 
it is natural to separate out the operators characterized by the minimum uncer- 
tainty relation ^(6 — /?)](& — 0)k) — $jk- Such operators are coherent vectors. 

Let us verify if vectors ipg of form (|1.4|) satisfy Equation. I|1.9I) in the Gaussian 
case. Substituting l|1.4|) . 1)2. 2f> into (|1.7|l and taking into account that 

(p\a) = exp {(3 f a - (a^a + f3 f (3)/2} 

we can find the conditional density 



p{fi\B) = J (p\a)(a\0)p(a\0)dii{a) 

= |JV + /| _1 e3tp{ - (/3-tf) t (iV + /)- 1 (/3-'!?)} ) (2.4) 

and also the function z(/3, i?) occurring in Equation. I|1.9|l : 

i(/3,d) =\n\l + S{N + iy^+^S-H - (-d - ApyG(d - A(3), (2.5) 

where A — S(S + N + I)- 1 , and G = S^ 1 + (N + J) -1 . In £3 only the 
last term is significant for Equation. 1)1. 9JI : the first two terms drop out after 
substitution into this equation. Substituting the last term (A/3— d)^ G(A{3— d) = 
c(A(3, i?) of function 1)2.5)1 into Equation. 1)1.9)1 and using representation 1)2.2)1 of 
the operator we verify that the coherent vector y>g = |/3) is a solution of 

Equation. ()1.9)) . Thus the solution of the formulated problem is reduced, as also 
in the optimization according to the mean square criterion, to the verification 
of the identity 

j y \a\a')(a'\(3) (c(A0,0) - J 

x p(a'\i9)dfi(a')P(di9) = 0, (2.6) 

where p(a|i?), P(d0) are defined in J221, 123 ■ 

Below we shall require the following formulas of integration in the complex 
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r-dimensional space C r 3 z: 

exp { - (z - a) t Q(z - (3)}\Q\ d^i(z) = 1, 
zexp{ - (z-a^Qiz- (3)}\Q\dfi(z) = (3, 

z* exp { - (z - a)^Q(z ~ /3)}\Q\ dfi(z) = a*, (2.7) 

(z - a)**H(z -/3)exp{ - (z - a)^ Q{z - [3)}\Q\d^{z) = SpQ^H, (2.8) 



all positive-definite (r x r)-matrices Q 7 H, where Sp is the trace in the space C 
where all vectors a,/3 £ C r are taking values. Taking into account that 

<a' |/3) </?'|/?)«a' I/?})- 1 = exp { - (/?' - a')\(3' - /?)}, 

m the identity J23 we can carry out the integration over after making the 
substitution A0 = z and making use of formula 1|2.7|) for Q = (AA^)^ 1 , H = G. 
As a result the expression in the parentheses in Ij2.6|) becomes {(3 — a'Y A^ G(A(3 — 
d)SpAA' r G. The integration of the last expression over a' with the density 

(a\a')(a'\[3}p(a'\ti) 

= (a\0)\N + I\~ l \Q\ exp {-(•&- a)*(N + I)" 1 ^ -(3)- {a' - a) + Q(a' - (3)} 

amounts to the use of the first and the third formulas of integration 12.71) . 
where for Q,a,/3 we should take Q = N^ 1 + I, a = (N + I) -1 (tf + Na), 
/3 = [N + I)- 1 ^ + N(3). This gives 

(P - (N + I)- 1 ^ + Na)) ] A^G{Af3 - 0) - Sp AA^l (a\0}\ 
x |JV + J|- 1 exp{ - (^-a) t (A^ + /)- 1 (i9-/3)}. (2.9) 

In order to complete the verification of identity ({2.6(1 it only remains to carry out 
averaging with the distribution P(di3). Rewriting the expression in the square 
brackets in l|2.9|) in the form 

(0 - (TV + I)(3 + NaYiN + J)~ X A^G{6 - A(3) - Sp AA^G 

and integrating it with respect to $ with the density 

\N + /rVr 1 exp { - (0 - a) f (iV + I)- 1 ^ -(3)- ^S- 1 ®} 
= \G\\S + N + I\e-^- Aa ^ G ^- A ^ 

according to the formulas (|2.7jl for Q = G,H = (N + iy 1 A^G we obtain 

\S + N + I| _1 Sp (G _1 (JV + I)~ l A^ - AA^)G. 



7 



This expression coincides with the left-hand side of 1)2. 6fl with an accuracy 
up to a nonzero factor (a\/3). Due to the identity 

G- 1 = (S^iN + I)- 1 )' 1 = S(S + N + J) _1 (iV + I) = A(N + I), 

as a consequence of which 

Sp (GT^JV + fT 1 ^ - = Sp (AA* - AA+)G = 0, 

the validity of the verified equation is follows straightforward. 

Thus in a linear Gaussian Bosonic channel the optimum decoding is given 
by coherent vectors of the form (|1.4|) . The ideal quasi- measurement, defined by 
the overcomplete nonorthogonal set {\(3)} of coherent vectors, is realized by the 
measurement of complex amplitudes 

fi = b + a+ = $ + a + a^ = d + a, (2.10) 

in the expanded space these amplitudes are of a complete orthonormal set of 
eigenvectors. The measurement of simultaneous observables H2.10J) is a linear 
indirect measurement of noncommuting b — (bj) investigated in £Q. The max- 
imum amount of information decoded by the optimal ideal quasi-measurement 
is easily found by averaging the function l|2.5|l : 

I = ln|/ + 5'(iV + /)- 1 | = Sp In (I +S{N + 1)- 1 ). 

In the one-dimensional case the obtained amount of information coincides with 
that shown in [SJ EI ■ 

In conclusion, we note that if a Gaussian stationary signal transmitted along 
the channel is mixed with stable thermal noise of temperature 8 = /cT, then 
the rate information transmission in optimum decoding is determined by the 
information 

In (l + s v /(n v + 1)) 
transmitted by each mode summed over the operating frequency range N 3 v. 

1= f ln(l + s v (n v + l)- 1 )dv = f {l + s v (l-e- hv l 9 ))dv. (2.11) 

JN JN 

Here s v is the spectral intensity of the transmitted signal expressed in units 
of hv, and n v = (e hv l B — 1) _1 is the spectral intensity of the noise. In the 
classical limit hv/6 <C 1 the rate of information transmission goes over into 
the corresponding classical expression I = Jln(l + <x„/ff)dv, where a 2 v — s v hv. 
In the opposite case hv/9 1 of low temperatures the transmission rate I = 
Jln(l + s u )Av remains finite in contrast to the classical case. It is not difficult 
to see that the quantum corrections to (|2.11|) are significant only for weak (at 
the receiving end of the channel) signals <C 1, for which formula (|2.11|l can 
be written in the form I » J(l — e~ hh '/ e )s u dv . 
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